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£ ■ ABSTRACT 

We have used the VLA at 43 GHz to image the radio continuum emission 
j> ■ from o Ceti, R Leo, and W Hya and to precisely locate their SiO maser emission 

with respect to the star. The radio continuum emission region for all three stars 
t— i ■ has a diameter close to 5.6 AU. These diameters are similar to those measured 

at infrared wavelengths in bands containing strong molecular opacity and about 
twice those measured in line-free regions of the infrared spectrum. Thus, the 
radio photosphere and the infrared molecular layer appear to be coextensive. 
The 43 GHz continuum emission is consistent with temperatures near 1600 K 
and opacity from H~ free-free interactions. While the continuum image of o Ceti 
appears nearly circular, both R Leo and W Hya display significant elongations. 
The SiO masers for all three stars show partial rings with diameters close to 
8 AU. 

Subject headings: Stars: Long-Period Variables, Circumstellar Matter, Atmo- 
spheres - Radio Continuum: Stars - Masers 
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1. Introduction 



A typical oxygen-rich Mira variable has a radius, i?* , defined by line-free regions of the 
spectrum, of about 1-2 AU. Surrounding the star is an extensive "molecular layer," which 
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extends to nearly 2R* and can be opaque across most of the visible and infrared spectrum 
The time-variable formation of metallic oxides in the molecular photosph ere can explain the 



thou sand-fold (8 mag) variations in the optical light curves of some Miras (IReid &: Goldston 



20021 ). Also at about 2/2* one detects a "radio photosphere," surrounded by an SiO maser 
emitting region. Beyond the SiO region dust grains form, radiation pressure accelerates the 
dust, transferring outward momentum to the gas and leading to significant mass loss. 



Reid fc Menten I (jl997al ) published measurements made with the Very Large Array 



(VLA) of a sample of Mira-like variables (o Ceti, R Leo, W Hya, R Aql, x Cyg, and R Cas). 
Their results revealed that the these stars exhibit a cm- wavelength spectral index of ~ 1.9 
and nearly constant radio light curves (variation ^ ±15% about the mean). Also, they 
imaged the radio continuum of W Hya at 22 GHz and found the star precisely at the center 
of a ring of H2O masers with a radius of 0.15 (IReid fc Menten Ill990l ). 



The detection of a "radio photosphere" near 2R* provides a new observational tool to 
increase our understanding of the region between the star and the dust formation zone. For 
the physical conditions expected in the radio photosphere, free electrons, obtained mostly 
from the ionization of potassium and s odium, provide the dom inant opacity through free- 
free interactions with neutral H and H2. IReid fc Menten I (Il997al ) proposed a simple physical 
model which can explain the radio frequency observations. At cm-wavelengths, unity optical 
depth is achieved at a radius of about 3 to 4 AU, where the density and temperature are 
~ 10 12 cm -3 and « 1600 K, respectively. This produces a flux density which is roughly 
twice that expected for black-body emission from the underlying star, as is observed. 

In this paper, we report observations at 7 mm wavelength (43 GHz) in which we simul- 
taneously imaged the radio photospheres and SiO maser emission from two Miras (o Ceti and 
R Leo) and one semi- regular variable (W Hya). These observations and supporting models 
constrain the physical and dynamical conditions between 2 and 3 R* 



2. Observations & Data Analysis 

The observations were made with the NRAO Very Large Array (VLA) under program 
AR446. Our goal was to measure directly the angular size of radio photospheres with the 
highest angular resolution afforded by the VLA. Our synthesized beam at 43 GHz in the 
A-configuration was ~ 0''04 full-width at half- maximum (FWHM), which well resolves a 



1 The National Radio Astronomy Observatory (NRAO) is operated by Associated Universities, Inc., under 
a cooperative agreement with the National Science Foundation. 
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typical Mira to a distance of ~ 200 pc. Since radio frequency "seeing" is limited primarily 
by fluctuations of water vapor in the atmosphere and is rarely bet ter than 0''1, we used a 



calibration scheme developed by us (IReid &: Menten I Il990l . Il997al ) in which circumstellar 



masers provide a phase reference for the continuum data. This allows diffraction limited 
imaging to be achieved. 

Between 2000 October 25 and November 6, we observed three stars: o Ceti, R Leo, and 
W Hya. Optical stellar phases for the three stars were approximately 0.05, 0.55, and 0.25, 
based on examination of AAVSO visual magnitudes for several cycles before and after the 
observations. 

The observations to measure stellar sizes were conducted in 4- IF continuum mode. The 
A-C IF pair (50 MHz bandwidth in right and left circular polarization) was centered at 
43164.9 MHz, which is 43 MHz above the rest frequency of 43122.08 MHz of the SiO v= 1, 
J=l-0 line. The B-D IF pair used a 1.563 MHz bandwidth centered on the SiO line at local 
standard-of-rest (LSR) velocities of 47.0 km s _1 for o Ceti, 1.0 km s" 1 for R Leo, and 
42.0 km s" 1 for W Hya. 

The interferometer data were calibrated using the Astronomical Image Processing Sys- 
tem (AIPS). Absolute flux density calibration was obtained fr om an observation of 3C286, 



assuming a flux density extrapolated fro m the formula given bvlBaars et al. I (119771 ). Follow- 



ing the procedures described in detail in IReid fc Menten I (ll997af ). we used roughly hourly 



observations of a strong continuum source to update pointing, to remove the electronic phase 
differences among the bands and to monitor the slowly-varying amplitude response of the an- 
tennas. The narrow-band data, containing very intense SiO maser emission, were phase-only 
self- calibrated on each 10 sec integration. The resulting antenna-dependent phase correc- 
tions were then applied to the broad-band data containing the continuum emission from the 
star. Since the masers are observed simultaneously with the continuum data and come from 
essentially the same position on the sky, this calibration removes all significant atmospheric 
phase-fluct uat ions . 

Once the data were "cross-self-calibrated," both the narrow-band maser and broad-band 
continuum signals were imaged with the AIPS task IMAGR. Map resolutions and noise levels 
are given in Table [TJ In order to simplify map analysis, CLEAN maps used a round restoring 
beam equal approximately to the geometric mean of the dirty beam's major and minor axes. 
These maps are shown in left column of panels in Fig. [TJ 

In addition to images, we also analyzed the interferometer visibility data in the following 
manner. Firstly, we measured the position of peak continuum emission by fitting an elliptical 
Gaussian brightness distribution to the map. Next, we "shifted" the data to place the 
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continuum peak at the interferometer phase-center. After time averaging the visibilities for 
1 minute, we fitted the data in the (u,v)-p\a,ne with a uniform-brightness, elliptical-disk 
model. This fitting was done using MIRIAD analysis software. Parameters of these fits are 
given in Table [2J 

Displaying the 2-dimensional (u, ■uj-plane data and model fits for weak sources is difficult. 
In order to obtain sufficient signal-to-noise ratio, the data need to be vector averaged over 
large portions of the irregularly sampled (u,v)-Y>\&ns. This precludes, for example, plotting 
visibility amplitude versus projected baseline length along a narrow strip at given position 
angle in the (u,v)-p\a.ne. We have chosen to bin the data and plot the real part of the 
visibility versus baseline length. These plots appear in the right-hand panels of Fig. [TJ For 
an elliptical source, this mixes data with differing visibilities, owing to the greater resolution 
of a source for a given baseline length if it is aligned with the source major axis than the 
minor axis. In order to indicate the magnitude of this effect, we have plotted the theoretical 
visibilities for circular sources with diameters equal to the major and minor axes of the fitted 
elliptical model. 

We also observed the SiO maser emission at high spectral resolution with several scans 
in spectral-line mode interspersed among the dual-band continuum observations. We covered 
the entire SiO line emission with a bandwidth of 3.125 MHz and used 64 spectral channels, 
which provided a velocity resolution of 0.34 km s _1 . The line data were self-calibrated by 
choosing a channel with strong emission as a reference, and the resulting phase and amplitude 
corrections were applied to the other channels. Because the maser emission was very strong 
(~ 10 3 Jy), we did not require bandpass calibration. From the line data we produced a 
spectral-line data cube, which we restored with a 20 mas FWHM Gaussian beam, taking 
advantage of the high signal-to-noise ratio and nearly point-like emission in each spectral 
channel to "super-resolve" by a factor of about two. 

The continuum emission (in the 50 MHz A-C IF pair) and maser emission (in the 1.563 
MHz B-D IF pair), which were observed simultaneiously, are automatically aligned by the 
calibration procedures to an accuracy of better than 1 mas. In order to align the maser 
emission from the spectral-line data, not observed simultaneously, with the continuum, we 
performed the following steps: 1) we produced a pseudo-continuum map from the line data, 
using spectral channels covering the velocities that were within the 1.563 MHz passband of 
the dual-band continuum setup; 2) we overlaid the pseduo-continuum map and on the map 
from the 1.563 MHz passband; and 3) we translated the pseduo-continuum map to achieve 
the best alignment by eye. Since the maps were nearly identical and contained unresolved 
spots with high signal-to-noise ratios, alignment of the continuum and maser emission to 
better than « 2 mas (or 10% of the restoring beam) was achieved. Based on this cross- 
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registration, we show velocity-integrated SiO maser emission superposed on the continuum 
emission of the stars in Fig. [2j In order to avoid introducing excess noise in the velocity- 
integrated SiO maser maps, we set the flux density of all pixels with less than 1 Jy beam -1 
to zero before summing the emission of the individual spectral channels. 



3. Radio Photospheric Characteristics 



3.1. Diameters 



We clearly resolve the radio photospheres of all three stars. For o Ceti, we find a 
nearly circular image with a uniform-disk diameter of about 52 mas (see Table El) . At a 
distance of 110 pc, based on the Period-Luminosity (PL) relation of iFeast et al. I (119891 ) and 
the photometry of IHaniff. Scholz fc Tuthill I (119951 ). this corresponds to a radio diameter 
of 5.7 AU. Throughout this paper, we adopt the PL-based distances, which appear to be 
more accurate than Hipp arcos parallaxes for these large, variable and irregularly bright stars 
(IReid fc Menten Ill997al ). We measure a mean diameter of 50 mas for R Leo, corresponding 
to 5.5 AU at a distance of 110 pc. W Hya displays a mean angular diameter of 58 mas, 
which at 95 pc distance implies a physical diameter of 5.5 AU. Thus, the radio-frequency 
mean diameters of all three stars are remarkably similar, averaging close to 5.6 AU. 

The diameters of the radio photospheres are similar to those measured in the visible and 
infrared in narrow bands that contain strong molecular opacity. For example, in a 10-Atm- 



wide b and centered at 710 fim containing strong TiO absorption, IHaniff. Scholz &: Tuthill 



( 119951 ) find (uniform-disk) ang ular diameter s of 53 , 45 and 66 mas for o Ceti, R Leo and 
W Hya, respectively. Similarly, iPerrin et al. I (120041 ) infer a "molecular layer" diameter of 52 
and 50 mas for o Ceti and R Leo, respectively. Generally it is difficult to "see" through the 
molecular photosphere and down to R* . Using n arrow bands placed in spectral regions with 
minimal molecular opacity, IPerrin et al. I (120041 ) find stellar diameters of 26 mas for o Ceti 
and 22 mas for R Leo. Model-fit Rosseland mea n diameters from observations using broade r 
bands yield diameters that are ~ 20% larger (jWoodruff et al. I l2004t iFedele et al. 1 120051 ) . 
Thus, the radio and molecular layers occur at radii of « 2/?*, adopting the narrow-band 
stellar diameters. 



3.2. Mira Binary 



Mira A is o Ceti, an M-type giant, and Mira B is possibly a low-mass main-sequence 
star ([Ireland et al. 1 120071 ) or a white dwarf. [Matthews Karovska I (120061 ) detected both 
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stars in the Mira binary with the VLA. They found the 43 GHz flux density of Mira B is 
weaker than that of o Ceti by a factor of ~ 3. We well resolve the binary, which has a 
projected separation of ~ 500 mas, and fail to detect Mira B. Our 3a limits are 0.4 and 
0.7 mJy beam -1 for untapered and tapered maps with (geometric mean) interferometer 
beams of 50 and 250 mas, respectively. Were the radio emission from Mira B unresolved and 
at the ~ 0.9 mJy level seen by Matthews & Karovska, we would have detected it. Since the 
Matthews & Karovska B-configuration observations at 43 GHz had a beam size comparable 
to our tapered map, it appears that angular resolution alone cannot explain our failure to 
detect Mira B and that it must be time variable. Since we well resolve (~ 10 beams) the 
binary and Mira B is not detected, our o Ceti results should not be significantly affected by 
the secondary. 



3.3. Asymmetries 

While o Ceti's image is circular to within ±4%, we detect non-circular radio photo- 
spheres for R Leo and W Hya. The differences between the major and minor axes are 
significant at about the 2a level. We are reasonably confident that the non-circular struc- 
tures do not come from instrumental limitations for several reasons: 1) for R Leo the intrinsic 
interferometer beam is nearly round and for W Hya the image is elongated nearly perpen- 
dicular to the intrinsic interferometer beam; 2) analysis of the data both in the image and 
(u,v)-p\&ne yield similar results; and 3) the plots of visibility versus (u,v)-distance for the el- 
liptically shaped stars show scatter larger than expected for a circular source, but consistent 
with an elliptical source. 



3.4. Physical Conditions 

Both o Ceti and R Leo have (uniform-disk) brightness tempera tures near 1650 K. Thes e 



temperatures are close to those expected for the radio photospheres (IReid fc Menten Ill997al ) . 
We measure a somewhat higher brightness temperature of 2380 ±550 K for W Hya. Since the 
collisional ionization of metals, and hence the H~ free-free opacity, is strongly temperature 
dependent below 2500 K, the higher measured brightness temperature of W Hya might 
suggest that the total density in the W Hya photosphere is lower than in o Ceti or R Leo. 
However, the uncertainties in brightness temperature are fairly large, owing predominantly 
to the squared-dependence on the uncertainty in the size measurement, and all three stars 
could have similar temperatures. Our temperature measurements are consistent, within 
uncertainties, with our earlier measurement of brightness temperature of W Hya of 1570 ± 
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570 K at a lower frequency of 22 GHz Jlleid fc Menten Ill997af l. 



The brightness tempe ratures of the radio photospheres of o Ceti and R Leo are also 
similar to those inferred by lPerrin et al. I (120041 ) of about 2000 and 1600 K for the molecular 
layers of these stars. Spectroscopic observations of H2O in the infrared also indicate temper 
atures and densities close to those we previously measured and modeled. iMatsuura et al. 
(120021 ) model their 2.5 to 4.0 /im ISO telescope data with a hot (2000 K) layer close to the 
stellar photosphere (between 1 and 2R* ) and a cool (1400 to 1000 K) layer between 2.5 and 
4.R* . The total hydrogen densities in these layers are estimated to be ~ 10 12-13 cm~ 3 . These 
layers roughly straddle the radio photosphere where we hav e previously found temp eratures 



of about 1600 K and hydrogen densities of ~ 10 cm (IReid fc Menten Ill997al ). Since 



both the radii, temperatures and densities measured in the radio and infrared agree, this 
provides added confidence in both measurements. 



SiO Masers 



The SiO masers in o Ceti, R Leo and W Hya form partial rings with radii of about 40, 35 
and 41 mas, respectively (Fig. [2]). These radii were measured by visually overlaying circles 
of varying radii centered on the stellar continuum position and are probably accurate to 
±10%, limited mostly by the lack of complete maser rings. These partial rings are consistent 
with, but per h aps s lightly larger than, the Very Long Baseline Array (VLBA) results by 
Cotton et al. I (120041 ). who find radii at different epochs ranging from 30-38, 26-30, and 
42 mas for o Ceti, R Leo and W Hya, respectively, for the SiO v= 1, J=l-0 transition . The 
differences in ring radii between our VLA images, which are sensitive to lower brightness 
emission, and the VLBA results, which have higher angular resolution and use a d i fferen t 
approach to define the ring radii, are small. We note that iDiamond fc Kemball 63), 



in their extensive VLBA monitoring of TX Cam, find lower brightness SiO masers radially 
outside of the brightest masers. 

When we convert our ring diameters from angular to linear units, we find radii of 4.4, 3.8 
and 3.9 AU, respectively. Thus, these three stars show very similar SiO masing radii close to 
AU. Previous observations have also detected partial rings of SiO masers in Mira-like stars 



(|Mivoshi et al. 1 11994 iGreenhill et al. Ill995t iBoboltz & Marvel Il2000l : Diamond & Kemball 



20031 ) and have usually assumed that the stars were near the center of the rings. However, 
for the first time, we are able to demonstrate that the maser rings are precisely centered 
on the stars. These results nicely complement our previous result that the H2O masers 
in W Hy a also exhibit a partia l ring and are precisely centered on the stellar continuum 
emission (IReid & Menten Ill990l ). We plot the integrated H 2 maser emission from W Hya 
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in Fig. [31 superposed on its SiO maser and radio continuum emission. The H2O maser ring 
had a radius of ~ 150 mas, which is about 3.7-times that of the SiO masers. 

The simplest interpretation for the SiO masers is that they form in a spherical shell 
near radii of about 3-R* . If one adopts the somewhat smaller SiO radii from VLBA maps, 
which sample the higher brightness masers, and the larger stellar radii, from Rosseland mean 
opacities, the maser shells would have radii closer to 2R* . The inner boundary of the SiO 
masers is det ermined either by densiti es 10 10 cm -3 , leading to thermalization of the level 



populations (ILockett fc Elitzur 



1992 



f), o r by high radio-continuum opacity from H free- 



free interactions (IReid fc Menten Ill997al ). The outer boundary occurs when SiO densities 
become too low for strong maser amplification. Note that in addition to the rapid decrease 
in total density with increasing radius in the extended atmosphere, the SiO density can fall 
even faster owing to silicate grain formation. 



Proper motions of SiO masers in TX Cam by iDiamond fc Kemball I (120031 ) indicate 
complex dynamics. Expansion at a speed of ~ 7 (-D/390 pc) km s -1 , where D is the 
star's distance, is the dominant motion. However, some SiO maser spots are also observed 



to move inwards at a 



ew km s 1 , suggesting a possible large-scale convective structure. 



Diamond &: Kemball I (120031) adopt a distance ba s ed on the period-luminosity and photom- 



etry compiled by lOlivier. Whitelock fc Marang I (120011 ). An alternative distance estimate 
could come from the radius of the SiO maser ring, since we find all three of our stars have a 
ring radius near 4 AU. Adopting this ring radius for TX Cam, the measured mean angular 
radius of 16 mas gives a distance of 250 pc. At this distance the dominant SiO expan- 
sion speed would only be 4.5 km s _1 . In any event, there are no indications of pervasive, 
fast-moving 10 km s _1 ) shocks in the SiO masing region. 



Discussion 



A simple, consistent picture of the inner envelopes of oxygen-rich Mira variables is 
emerging from radio and infrared observations. Radio observations of stars like o Ceti, R Leo 
and W Hya provide data which constrains the kinematics and physical conditions between 
2 and 3i?*. Radio continuum emission comes predominantly from H~ free- free interactions. 
An optical depth of unity is observed at ~ 2R* and requires a te mperature and density of 



1600 K and ~ 10 12 cm" 3 , respectively (IReid k Menten Ill997af ). 



The radio photosphere corresponds very closely to the molecular layer, identified at 
i nfrared wavelengths. The molecular layer can be directly observed in lin es of H2O and CO 
Jrlaniff. Scholz fc Tuthilllll995l : iMatsuura et al. II2OO2I ; IPerrin et aul2004h . As the star cools 
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cyclically, this region can become nearly opaque at visible wavelengths, owing to metallic 
oxide formation. This leads to dramatic changes in visible light as the temperature in the 
molecular photosphere at « 2R* can be low enough for visible light to "fall off" the Wien 
side of the black body spectrum (IReid &: Goldston II2002I ). 



SiO maser emission surrounds the radio/molecular layer and appears as partial rings 
with radii of 2 to 3.B*. These rings sometimes disp lay significant ellipticity as, for example, 
shown in TX Cam by iDiamond fc Kemball I (120031 ). SiO masers trace both outwardly and 
inwardly moving material, possibly indicative of large scale convection. Temperatures and 
densities inferred from SiO masers are ~ 1000 K and ^ 10 10 cm -3 . SiO maser proper motions 
suggest that outside of 3 R* shocks are mostly damped and propagate at ^ 5— 7 km s _1 . This 
complements a similar conclusion, that strong, fast shocks are not found i n the underlying 



radio ph otosphere, based on the observed nearly constant radio light curves (IReid fc Menten 



1997a| jbl). Thus, while radio data indicate a dynamic region between 2 and 3 R*, there are 
no signs of strong shocks heating large portions of the atmosphere to ~ 10, 000 K near 2R* 
as in some models of these stars (IWillson Il2000l ). 



The "radio photospheres" of o Ceti and W Hya appear somewhat elongated. The non- 
circular structure we detect in the radio photospheres of two AGB stars may be of significance 
in understanding the highly elongated structures seen in planetary nebulae, following a post- 
AGB phase. Clearly second epoch images of the radio photospheres, after several pulsation 
cycles, would be very interesting. If the non-circular structures change significantly, either in 
shape or orientation, this would suggest a (non-radial) pulsation origin. Alternatively, should 
the structures be stable over time, this would suggest a possible rotational or magnetic field 
origin, which may play a role in creating asymmetric structures in planetary nebulae. 

We thank L. Matthews for comments on an early version of the manuscript. We acknowl- 
edge with thanks the variable star observations from the AAVSO International Database 
contributed by observers worldwide and used in this research. 
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Table 1. Interferometer Parameters 



Star 


9 B (maj) 
(mas) 


0s(miri) 
(mas) 


P.A. 

(deg.) 


Continuum RMS Noise 
(mjy) 


Line RMS Noise 
(mjy) 


o Ceti 


53 


39 


-13 


0.13 


10 


R Leo 


61 


39 


-20 


0.14 


15 


W Hya 


81 


37 


- 7 


0.30 


20 



Note. — Beam major, 9g(maj), and minor, 9g(min), axes and East-of-North position angles, P.A., are from the best fit 
elliptical Gaussians to the "dirty beams." The maps displayed in Fig. [T] and Fig. |21 used round restoring beams as discussed in 
the text. 



Table 2. Measured Stellar Parameters 



Star 


Distance 


Opt. Phase 






P.A. 




T b 


SiO Radius 




(pc) 




(mas) 


(mas) 


(deg.) 


(mJy) 


(K) 


(mas) 


o Ceti 


110 ±9 


0.05 


54 ± 5 


50 ±5 


+39 ± 50 


4.8 ± 0.2 


1680 ± 250 


40 + 4 


R Leo 


110 ±9 


0.55 


61 ± 10 


39 ±6 


-20 ± 12 


4.1 ± 0.2 


1630 ± 410 


35 + 4 


W Hya 


95 ±8 


0.25 


69 ± 10 


46 ±7 


+83 ± 18 


8.0 + 0.4 


2380 ± 550 


41 + 6 



Note. — Distances are from lHaniff. Scholz k, Tuthill I l|l995h based on the PL relation of lFeast et al. ] l|l989h . Stellar major, 
Smaj, and minor, 9 m i„, diameters and East-of-North position angles, P.A., are from best-fit, uniform-brightness (Tj,), elliptical- 
disk models. Flux density, S„, uncertainties do not include a roughly ±10% uncertainty in the absolute flux density scale. SiO 
maser shell radii are estimated from the brightness peaks in Fig. [3] 
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Fig. 1. — Images shown as contours (left panels) and visibility vs. interferometer baseline 
length (right panels) of the 43 GHz continuum emission from o Ceti (top panels), R Leo 
(middle panels) and W Hya (bottom panels). Contour levels are integer multiples of 0.5 
mJy for o Ceti and R Leo and multiples of 1.0 mJy for W Hya. In the image panels, the 
interferometer FWHM "dirty" beams are shown as solid ellipses and the round restoring 
beams are shown as dotted circles in the lower left. The dashed ellipses superposed on 
the images are the sizes of uniformly bright elliptical disks that best fit the data. In the 
visibility plots, the two model curves correspond to expected visibilities for circular disks 
with diameters corresponding to the major and minor axes listed in Table EJ 




Fig. 2. — SiO maser emission (contours and gray scale) superposed on the fitted uniformly 
bright elliptical disk models of the radio continuum emission (green filled ellipses) from o Ceti 
(top panel), R Leo (middle panel) and W Hya (bottom panel). SiO contour levels are 1, 2, 4, 
8, 16, and 32 times 100, 40, and 200 Jy beam -1 for the three stars, respectively. The relative 
positions of the SiO masers and the radio continuum have been measured with accuracies of 
about 2 mas. The restoring beams for the radio continuum (larger dotted circle) and SiO 
maser data (smaller dotted circle) are indicated in the lower left corner of each panel. The 
true stellar sizes (re d filled circles) are from circular, uniformly bright, disk models of IR 
interferometer data (IPerrin et al. 1 12004 ) . 
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Fig. 3. — W Hya H 2 maser emission {gray scale) from iReid fc Menten I (11990l ). superposed 
on the SiO maser emission {yellow contours) and the radio continuum emission (elliptical 
disk model) {green filled ellipse) from Fig. [21 The H 2 velocity-integrated flux scale in 
Jy beam -1 is shown along the right side of the figure. The relative positions of the H 2 
masers and the radio continuum is accurate to about 10 mas. The restoring beam for the 
H 2 emission is indicated in the lower left corner. 



